Dynamics of glass phases in the two-dimensional gauge glass model 
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Large-scale simulations have been performed on the current-driven two-dimensional XY gauge 
glass model with resistively-shunted-junction dynamics. It is observed that the linear resistivity at 
low temperatures tends to zero, providing strong evidence of glass transition at finite temperature. 
Dynamic scaling analysis demonstrates that perfect collapses of current-voltage data can be achieved 
with the glass transition temperature Tg = 0.22, the correlation length critical exponent v — 1.8, 
and the dynamic critical exponent z — 2.0. A genuine continuous depinning transition is found 
at zero temperature. For creeping at low temperatures, critical exponents are evaluated and a 
non-Arrhenius creep motion is observed in the glass phase. 

PACS numbers: 74.25.-q, 68.35.Rh, 64.70.Q-, 05.10.-a 
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I. INTRODUCTION 

The evidences to support the existence of a vortex glass 
(VG) phase in strongly disordered type-II superconduc- 
tors have been reported in many experiments by the dy- 
namic scaling of the measured current- voltage data[l|. 
Theoretically, the XY gauge glass modeljl, y| is often 
used to describe the VG phase, although it lacks some 
of properties and symmetries due to the absence of net 
magnetic fields [1, 0, Q. Now there is general consen- 
sus that a finite-temperature VG transition occurs in the 
three-dimensional gauge glass model 

The situation is, however, much less clear in two di- 
mensions (2D). The experimental quest of the VG tran- 
sition in high- Tc cuprate films has provided contin- 
uous excitement and puzzles. Recently, in a positionally 
disordered Josephson junction arrays with the maximal 
disorder strength 10], where the 2D gauge glass model 
is realized, a possible finite-temperature glass transition 
has been observed experimentally. On the theoretical 
side, the existence of a finite-temperature glass transi- 
tion in the 213 gauge glass model remains a topic of 
controversy [II [Ij, Q [IE [11, [13, [11, [ll [3 ■ It is pre- 
dicted that, in the zero-temperature numerical renormal- 
ization group studies of domain walls and the calculations 
of stiffness exponents, there is no ordered phase at any 
finite temperature in 2D[i2, [l3|. On the other hand, the 
finite-temperature glass transition {Tg « 0.2J) has also 
been spported by extensive resistively-shunted-junction 
(RSJ) dynamic simulations [l6l . [TtI . [Hi and Monte Carlo 
simulations (20| . 

The depinning transition at zero temperature and 
the creep motion at low temperatures have attracted 
considerable attention both analytically [U [13, [13 and 
nume ricallylH, [H, Hi] in a large variety of physical sys- 
tems, such as charge density waves in solids, field-driven 
motion of domain walls in ferromagnets and flux lines in 
type-II superconductors. Since the non- linear dynamic 
response in these systems produces a rich physical pic- 
ture, there has been increasing interest in studies of these 
phenomena. 



In this paper, based on the RSJ dynamics, we perform 
large-scale dynamic simulations on the 2D gauge glass 
model. Both the glass transition temperature Tg and the 
critical exponents are estimated. The depinning tran- 
sition at zero temperature and the creep motion below 
Tg are also investigated. The rest of the paper is orga- 
nized as follows. Sec. II describes the model and dynamic 
method. Sec. Ill presents the main results, where some 
discussions are also made. Finally, a short summary is 
given in the last section. 



II. MODEL AND DYNAMIC METHOD 

The Hamiltonian of the 2D gauge glass model is given 
byOi] 
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where the sum is over all nearest neighbor pairs on a 2D 
square lattice, specifies the phase of the superconduct- 
ing order parameter on grain i, Jq denotes the strength 
of Josephson coupling between neighboring grains, and 
the quenched variable Aij is distributed uniformly in the 
interval [— 7r,7r). The present simulations are carried out 
with the system size L = 128 for all directions. 

The RSJ dynamics is incorporated in simulations, 
which can be described as 
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where Jext.z is the external current which vanishes ex- 
cept for the boundary sites. The rjij is the thermal noise 
current with zero mean and a correlator {riij{t)riij{t')) = 
2crkBT5{t — t'). In the following, the units are taken of 
2e = Jq = h = a = ks ~ I- 

In the present simulations, a uniform external current 
Lj: along the x direction is fed into the system. The fluc- 
tuating twist boundary condition ^3, [ll] is applied in 
both directions. The supercurrent between sites i and j 
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FIG. 1: Log-log plots of I — R curves at various temperatures. 

is now given by J-^j = Jq sm{9i — 9j — Aij — • A) 
with 9i — (pi + Ti ■ A. and A = (A^;, Aj^) the fluctuat- 
ing twist variable. The new phase angle 9i is periodic in 
both X and y directions. Then, the dynamics of A^ can 
be written as 

^"^^2 [Jt'-lj+Vi]]- Ia.a^x,y. (3) 

The voltage drop is F = —LK^. 

The above equations can be solved efficiently by a 
pseudo-spectral algorithm [29] due to the periodicity of 
the phase in all directions. The time stepping is done us- 
ing a second-order Runge-Kutta scheme with At — 0.05. 
The time-averaged voltages are calculated over a long- 
time scale after reaching the steady state. To deter- 
mine the steady state, we have checked (V)n for every 
(2" — 2"~^) time steps. We assume that the system 
reaches a steady state when the fluctuation of the mean 
voltage \{(y)n — (V)„_i)/(y)„| is less than 0.5% for sev- 
eral n's after n = 20. Once this criterion is satisfied, 
we record the (y)n as the final estimate of the voltage 
V . The value of n is typically 25 in the present simu- 
lations. The detailed procedure in the simulations was 
described in Ref. [2§|. We have performed simulations 
with 10 different realizations of disorder and observed 
that the results are quite close from sample to sample, so 
good self-averaging effects exist in the present large sys- 
tems. This point is also supported by a recent study in 
Josephson-junction arrays by Um et al. [l^. Our results 
below are averaged over 10 realizations of disorder. For 
the results presented in the following figures, error bars 
are smaller than or comparable with the symbol size. 

III. SIMULATION RESULTS AND 
DISCUSSIONS 

First, we study the possible glass phase transition. 
The glass transition temperature was estimated to be 
Tg w 0.20 by several groups [H [H, Hi]. So the current- 
voltage characteristics are measured at various tempera- 
tures ranging from 0.1 to 0.5, which cover the previous 
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FIG. 2: Dynamic scaling of I — R data at various temperatures 
according to Eq. (|4]). 

Tg. At each temperature, we try to probe the system at 
currents as low as possible. FiglT] displays the resistivity 
R = V/ 1 as a function of current / at various temper- 
atures in a log-log scale. It is obvious that, at lower 
temperatures, R tends to zero as the current decreases, 
suggesting that there is a true superconducting phase 
with zero linear resistivity. While at higher tempera- 
tures, R tends to a finite value, corresponding to Ohmic 
resistivity in the vortex liquid. These observations reveal 
strong evidence of the existence of the low-temperature 
glass phase in the 2D gauge glass model. 

Assuming that the VG transition is continuous and 
characterized by the divergence of the characteristic 
length and time scales i ~ ( z is the dynamic ex- 
ponent), Fisher, Fisher, and Huse [sO] proposed the fol- 
lowing dynamic scaling ansatz to analyze the glass tran- 
sition from a vortex liquid with Ohmic resistivity to a 
superconducting glass state 

j,j^^z+2-d ^ ^-^(/^'i-i/r), (4) 

where d is the dimension of the system (d = 2 in this 
paper), ^ oc| T/Tg — 1 j^'^ is the correlation length which 
diverges at the transition and ^'±(x) are scaling functions 
for T > Tg and T <Tg, respectively. Eq.© is often used 
to scale the measured current-voltage data in the VG 
transitions in experiments. 

To extract the critical behavior from the numerical re- 
sults of the current- voltage characteristics, we also per- 
form a dynamic scaling analysis. As shown in Fig [21 with 
Tg = 0.22 ±0.02 , z = 2.0 ±0.1, and v 1.8 ±0.1, an ex- 
cellent collapse is achieved according to Eq. ([4]) except for 
the curve of T = 0.1. The errors are estimated by tun- 
ing these critical values until the collapses become poor 
evidently. The curve at T = 0.1 is obviously beyond the 
critical regime. 

The finite-size effects are particularly significant at 
temperatures sufficiently close to Tg when the correla- 
tion length exceeds the system size. For the temperatures 
considered here and the very large system size L = 128, 
we believe that the finite-size effects are negligible in the 
present simulations. To confirm this point, we perform 
particular simulations right at Tg = 0.22 obtained above 
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FIG. 3: Dynamic scaling oi I — R data at Tg — 0.22 according 
to Eq.©. 



FIG. 4: Log-log plot of V versus (/ — Ic) curve at zero tem- 
perature. 



for different system sizes. At Tg, the correlation length 
is cut off by the system size in any finite system, so the 
scaling form Eq.Q becomes 



(5) 



A good collapse is illustrated in Fig[3]with z — 2.0. This 
consistency demonstrates that the results estimated from 
FiglH are reliable. Therefore new evidence of a finite- 
temperature glass transition is provided convincingly in 
the 2D gauge glass model. 

The obtained Tg and dynamic exponent z are well con- 
sistent with those in equilibrium RSJ simulations [l^ and 
Monte Carlo simulations [2^. The value of estimated 
here is larger than those ( 1.1 ^ 1.2) obtained in several 
simulations [l 7. 20], but still falls in the range of 1.0 ~ 2.0 
usually observed at the VG transitions experimentally. 
Within the same RSJ dynamics, the finite-size scaling 
for the linear resistivity for sample sizes L < 10 gives 
v = 1.2(2) in Ref.[l3]- The discrepancy may originate 
from the different scaling method and sizes used. In a 
previous conference paper [31 by one of the present au- 
thor and collaborators, with the uses of a different simu- 
lation approach and a scaling form different from Eq. ^ 
slightly by removing temperature T, Tg = 0.22 , z = 2.0, 
and v = 1.2 were obtained. The present simulation gives 
larger value of v. 

Based on the analysis of data at teniperatures above 
Tg — 0.22, previous RSJ simulations [isj with an open 
boundary condition in the 2D gauge glass model demon- 
strated a zero-temperature criticality. It has been shown 
that the voltage drop next to the boundary regime is 
particularly large and dominates the total voltage drop 
across the sample at low currents HH, 113 • There- 
fore, one should measure the voltage drop inside the 
sample [28l [s^ . So the conclusion based on the total 
voltage drop across the system with the open boundary 
condition[3, [l^ may not be reliable. 

Interestingly, in experiments on Nb wire networks [ssj. 
the critical exponents v = 1.7 ^ 1.9 were obtained for 
high filling factors / = 1/2,0.618, and 2/5, which are 
very close to the present value. It was suggested in Ref. 
that the superconducting state and transitions in 



the networks become independent of / in the gauge glass 
limit. The further work is needed to clarify the relation 
between the experimental observations and the present 
simulations. 

To shed some light on the nature of this low- 
temperature glass phase, we will study the depinning and 
creep phenomena. 

At zero temperature, we start from high currents with 
random initial phase configurations. The currents are 
then lowered step by step. The steady-state phase con- 
figurations obtained at higher currents are chosen to be 
the initial phase configurations of lower currents in the 
next step. It becomes more and more difficult to mea- 
sure the voltage with decreasing currents. In the vicin- 
ity of the critical current, a huge amount of computer 
time is consumed to get accurate results. Fig H] presents 
the current-voltage characteristics at T = in a log- 
log scale. We observe a continuous depinning transition 
with a unique depinning current's^, which can be de- 
scribed asV (X {I ~ Icf with Ic = 0.2165 ± 0.0005 and 
(3 — 1.892 ± 0.003. Note that the depinning exponent (3 
is greater than 1, consistent with the mean field studies 
of charge density wave models [33|. 

At low temperatures, the current-voltage characteris- 
tics are rounded near the zero-temperature critical cur- 
rent due to thermal fluctuations. An obvious crossover 
between the depinning and creep motion can be observed 
around Ic at lower accessible temperatures. In order to 
address the thermal rounding of the depinning transi- 
tion. Fisher [s^l first suggested to map this system to fer- 
romagnets in fields where the second-order phase tran- 
sitions occur. This mapping was later extended to the 
random-field Ising modelf24'| and flux lines in type-II 
superconductors [25]. If the voltage is identified as the 
order parameter, the current and temperature are equiva- 
lent to the inverse temperature and the field in ferromag- 
netic systems, respectively, analogous to the second-order 
phase transitions, a scaling relation among the voltage, 
current and temperature in the present model should fol- 
low the form 



T/(T, /) = T^/»S[T-^'^\l - Ic/I)], 



(6) 
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y = 0.0417exp(1.17a;), which is also plotted in the Fig[6] 
with a solid line. Note that the product of the two ex- 
ponents (58 describes the temperature dependence of the 
creep law. Interestingly, (35 ~ 1.81 deviates from unity, 
demonstrating that the creep law is a non-Arrhcnius 
type. At higher temperatures, say T > 0.1, deviations 
from the scaling relation are also observed in FigjHl which 
can be attributed to strong thermal fluctuations. The 
non-Arrhenius type creep phenomena only take place at 
low temperatures. 



FIG. 5: 




FIG. 6: Scaling plot of / 
according to Eq.®. 
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where S{x) is a scaling function with S{x 0) =const. 

It is implied in Eq.® that right at / — /c the volt- 
age shows a power-law behavior V{T,I — Ic) oc T^/*, 
providing a tool to determine the critical exponent 1/6. 
The log-log V — T curves are plotted in Fig. [5] at three 
currents around Ic- We can see that the critical current 
is between 0.21 and 0.22. The values of voltage at other 
currents within (0.21, 0.22) can be evaluated by quadratic 
interpolation. The square deviations from the power law 
can be calculated. The current at which the square de- 
viation is minimum can be considered as the critical cur- 
rent Ic — 0.2165 ± 0.0005, consistent with that obtained 
at zero temperature. The temperature dependence of 
voltage at the critical current is also exhibited in Figl5l 
yielding 1/6 = 1.046 ± 0.002. 

With the values of (3, 6 and Ic obtained above, ac- 
cording to the scaling relation Eq.®, a scaling plot 
of the simulated current-voltage data in a wide range 
of temperature is presented in FiglS] without any ad- 
justable parameter. A perfect collapse of the data for 
temperatures T < 0.10, far below Tg = 0.22, to a single 
curve for currents less than Ic is clearly shown. This 
collapse can be fitted well to an exponential function 



IV. SUMMARY 

We have performed large-scale dynamic simulations of 
the 2D gauge glass model within the RSJ dynamics. The 
strong evidence of the low-temperature glass phase is pro- 
vided in the dynamic sense. By the dynamic scaling anal- 
ysis, two perfect collapses of simulated current-voltage 
data are achieved with Tg = 0.22 ± 0.02 , z = 2.0 ± 0.1, 
and v = 1.8 ±0.1. The values of Tg and z are in agreement 
with those in the previous equilibrium Monte Carlo sim- 
ulations. While the value of i' is larger than that in liter- 
ature, which is, however, closer to that in experiments in 
the gauge glass limit. We have also studied the dcpinning 
transition at zero temperature and creep motion at low 
temperatures in detail. A genuine continuous depinning 
transition is observed at zero temperature. With the no- 
tion of scaling and the critical exponents obtained from 
the simulations at zero temperature and at the critical 
current, a perfect collapse of the current- voltage data at 
low temperatures is exhibited. The value of (36 deviates 
from unity and the scaling curve is fitted well by an expo- 
nential function, suggesting a non-Arrhcnius type creep 
motion in the glass phase of the 2D gauge glass model. 

It is worthy to note that in this model the current- 
voltage characteristics in the whole temperature regime 
below Tg can almost be described in the framework of two 
critical phenomena. One is the thermal rounding of the 
depinning transition, which is a second-order-like phase 
transition. The other is the celebrated VG transition. 
Further experimental and theoretical studies are needed 
to clarify the relation between the present observations 
and experimental findings. 
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